Long-range signaling by phosphoprotein waves arising from bistability in protein kinase cascades

A hallmark of protein kinase/phosphatase cascades, including mitogen-activated protein kinase (MAPK) pathways, is the spatial separation of their components within cells. The top-level kinase, MAP3K, is phosphorylated at the cell membrane, and cytoplasmic kinases at sequential downstream levels (MAP2K and MAPK) spread the signal to distant targets. Given measured protein diffusivity and phosphatase activities, signal propagation by diffusion would result in a steep decline of MAP2K activity and low bisphosphorylated MAPK (ppMAPK) levels near the nucleus, especially in large cells, such as oocytes. Here, we show that bistability in a two-site MAPK (de)phosphorylation cycle generates a novel type of phosphoprotein wave that propagates from the surface deep into the cell interior. Positive feedback from ppMAPK to cytoplasmic MAP2K enhances the propagation span of the ppMAPK wave, making it possible to convey phosphorylation signals over exceedingly long distances. The finding of phosphorylation waves traveling with constant amplitude and high velocity may solve a long-standing enigma of survival signaling in developing neurons

Bridging Time Scales in Cellular Decision Making with a Stochastic Bistable Switch

Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms. Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.Comment: 14 pages, 4 figure

Protein sequestration generates a flexible ultrasensitive response in a genetic network

Ultrasensitive responses are crucial for cellular regulation. Protein sequestration, where an active protein is bound in an inactive complex by an inhibitor, can potentially generate ultrasensitivity. Here, in a synthetic genetic circuit in budding yeast, we show that sequestration of a basic leucine zipper transcription factor by a dominant-negative inhibitor converts a graded transcriptional response into a sharply ultrasensitive response, with apparent Hill coefficients up to 12. A simple quantitative model for this genetic network shows that both the threshold and the degree of ultrasensitivity depend upon the abundance of the inhibitor, exactly as we observed experimentally. The abundance of the inhibitor can be altered by simple mutation; thus, ultrasensitive responses mediated by protein sequestration are easily tuneable. Gene duplication of regulatory homodimers and loss-of-function mutations can create dominant negatives that sequester and inactivate the original regulator. The generation of flexible ultrasensitive responses is an unappreciated adaptive advantage that could explain the frequent evolutionary emergence of dominant negatives

A novel, rapid method to compare the therapeutic windows of oral anticoagulants using the Hill coefficient

A central challenge in designing and administering effective anticoagulants is achieving the proper therapeutic window and dosage for each patient. The Hill coefficient, nH, which measures the steepness of a dose-response relationship, may be a useful gauge of this therapeutic window. We sought to measure the Hill coefficient of available anticoagulants to gain insight into their therapeutic windows. We used a simple fluorometric in vitro assay to determine clotting activity in platelet poor plasma after exposure to various concentrations of anticoagulants. The Hill coefficient for argatroban was the lowest, at 1.7±0.2 (95% confidence interval, CI), and the Hill coefficient for fondaparinux was the highest, at 4.5±1.3 (95% CI). Thus, doubling the dose of fondaparinux from its IC50 would decrease coagulation activity by nearly a half, whereas doubling the dose of argatroban from its IC50 would decrease coagulation activity by merely one quarter. These results show a significant variation among the Hill coefficients, suggesting a similar variation in therapeutic windows among anticoagulants in our assay

The edges of understanding

A culture's icons are a window onto its soul. Few would disagree that, in the culture of molecular biology that dominated much of the life sciences for the last third of the 20th century, the dominant icon was the double helix. In the present, post-modern, 'systems biology' era, however, it is, arguably, the hairball

Unlimited multistability in multisite phosphorylation systems

Reversible phosphorylation on serine, threonine and tyrosine is the most widely studied posttranslational modification of proteins (1, 2). The number of phosphorylated sites on a protein (n) shows a significant increase from prokaryotes, with n less than or equal to 7 sites, to eukaryotes, with examples having n greater than or equal to 150 sites (3). Multisite phosphorylation has many roles (4, 5) and site conservation indicates that increasing numbers of sites cannot be due merely to promiscuous phosphorylation. A substrate with n sites has an exponential number (2^n) of phospho-forms and individual phospho-forms may have distinct biological effects (6, 7). The distribution of these phospho-forms and how this distribution is regulated have remained unknown. Here we show that, when kinase and phosphatase act in opposition on a multisite substrate, the system can exhibit distinct stable phospho-form distributions at steady state and that the maximum number of such distributions increases with n. Whereas some stable distributions are focused on a single phospho-form, others are more diffuse, giving the phospho-proteome the potential to behave as a fluid regulatory network able to encode information and flexibly respond to varying demands. Such plasticity may underlie complex information processing in eukaryotic cells (8) and suggests a functional advantage in having many sites. Our results follow from the unusual geometry of the steady-state phospho-form concentrations, which we show to constitute a rational algebraic curve, irrespective of n. We thereby reduce the complexity of calculating steady states from simulating 3 times 2^n differential equations to solving two algebraic equations, while treating parameters symbolically. We anticipate that these methods can be extended to systems with multiple substrates and multiple enzymes catalysing different modifications, as found in posttranslational modification 'codes' (9) such as the histone code (10, 11). Whereas simulations struggle with exponentially increasing molecular complexity, mathematical methods of the kind developed here can provide a new language in which to articulate the principles of cellular information processing (12)

Monomeric Bistability and the Role of Autoloops in Gene Regulation

Genetic toggle switches are widespread in gene regulatory networks (GRN). Bistability, namely the ability to choose among two different stable states, is an essential feature of switching and memory devices. Cells have many regulatory circuits able to provide bistability that endow a cell with efficient and reliable switching between different physiological modes of operation. It is often assumed that negative feedbacks with cooperative binding (i.e. the formation of dimers or multimers) are a prerequisite for bistability. Here we analyze the relation between bistability in GRN under monomeric regulation and the role of autoloops under a deterministic setting. Using a simple geometric argument, we show analytically that bistability can also emerge without multimeric regulation, provided that at least one regulatory autoloop is present

Computational modelling of cancerous mutations in the EGFR/ERK signalling pathway

This article has been made available through the Brunel Open Access Publishing Fund - Copyright @ 2009 Orton et al.BACKGROUND: The Epidermal Growth Factor Receptor (EGFR) activated Extracellular-signal Regulated Kinase (ERK) pathway is a critical cell signalling pathway that relays the signal for a cell to proliferate from the plasma membrane to the nucleus. Deregulation of the EGFR/ERK pathway due to alterations affecting the expression or function of a number of pathway components has long been associated with numerous forms of cancer. Under normal conditions, Epidermal Growth Factor (EGF) stimulates a rapid but transient activation of ERK as the signal is rapidly shutdown. Whereas, under cancerous mutation conditions the ERK signal cannot be shutdown and is sustained resulting in the constitutive activation of ERK and continual cell proliferation. In this study, we have used computational modelling techniques to investigate what effects various cancerous alterations have on the signalling flow through the ERK pathway. RESULTS: We have generated a new model of the EGFR activated ERK pathway, which was verified by our own experimental data. We then altered our model to represent various cancerous situations such as Ras, B-Raf and EGFR mutations, as well as EGFR overexpression. Analysis of the models showed that different cancerous situations resulted in different signalling patterns through the ERK pathway, especially when compared to the normal EGF signal pattern. Our model predicts that cancerous EGFR mutation and overexpression signals almost exclusively via the Rap1 pathway, predicting that this pathway is the best target for drugs. Furthermore, our model also highlights the importance of receptor degradation in normal and cancerous EGFR signalling, and suggests that receptor degradation is a key difference between the signalling from the EGF and Nerve Growth Factor (NGF) receptors. CONCLUSION: Our results suggest that different routes to ERK activation are being utilised in different cancerous situations which therefore has interesting implications for drug selection strategies. We also conducted a comparison of the critical differences between signalling from different growth factor receptors (namely EGFR, mutated EGFR, NGF, and Insulin) with our results suggesting the difference between the systems are large scale and can be attributed to the presence/absence of entire pathways rather than subtle difference in individual rate constants between the systems.This work was funded by the Department of Trade and Industry (DTI), under their Bioscience Beacon project programme. AG was funded by an industrial PhD studentship from Scottish Enterprise and Cyclacel

Oscillatory stimuli differentiate adapting circuit topologies

This is the author accepted manuscript. The final version is available from Springer Nature via the DOI in this record.Biology emerges from interactions between molecules, which are challenging to elucidate with current techniques. An orthogonal approach is to probe for 'response signatures' that identify specific circuit motifs. For example, bistability, hysteresis, or irreversibility are used to detect positive feedback loops. For adapting systems, such signatures are not known. Only two circuit motifs generate adaptation: negative feedback loops (NFLs) and incoherent feed-forward loops (IFFLs). On the basis of computational testing and mathematical proofs, we propose differential signatures: in response to oscillatory stimulation, NFLs but not IFFLs show refractory-period stabilization (robustness to changes in stimulus duration) or period skipping. Applying this approach to yeast, we identified the circuit dominating cell cycle timing. In Caenorhabditis elegans AWA neurons, which are crucial for chemotaxis, we uncovered a Ca2+ NFL leading to adaptation that would be difficult to find by other means. These response signatures allow direct access to the outlines of the wiring diagrams of adapting systems.The work was supported by US National Institutes of Health grant 5RO1-GM078153-07 (F.R.C.), NRSA Training Grant CA009673-36A1 (S.J.R.), a Merck Postdoctoral Fellowship at The Rockefeller University (S.J.R.), and the Simons Foundation (S.J.R.). J.L. was supported by a fellowship from the Boehringer Ingelheim Fonds. E.D.S. was partially supported by the US Office of Naval Research (ONR N00014-13-1-0074) and the US Air Force Office of Scientific Research (AFOSR FA9550-14-1-0060)

Phenotypic Variation and Bistable Switching in Bacteria

Microbial research generally focuses on clonal populations. However, bacterial cells with identical genotypes frequently display different phenotypes under identical conditions. This microbial cell individuality is receiving increasing attention in the literature because of its impact on cellular differentiation, survival under selective conditions, and the interaction of pathogens with their hosts. It is becoming clear that stochasticity in gene expression in conjunction with the architecture of the gene network that underlies the cellular processes can generate phenotypic variation. An important regulatory mechanism is the so-called positive feedback, in which a system reinforces its own response, for instance by stimulating the production of an activator. Bistability is an interesting and relevant phenomenon, in which two distinct subpopulations of cells showing discrete levels of gene expression coexist in a single culture. In this chapter, we address techniques and approaches used to establish phenotypic variation, and relate three well-characterized examples of bistability to the molecular mechanisms that govern these processes, with a focus on positive feedback.